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ABSTRACT 

We present a new high-mass membership of the nearby Sco 0B2 association based 
on HIPPARCOS positions, proper motions and parallaxes and radial velocities taken 



from the Kliarchenko et al. (12007) catalogue. The Bayesian membership selection 
method developed makes no distinction between subgroups of Sco 0B2 and utilises 
linear models in calculation of membership probabilities. We select 436 members, 88 
of which are new members not included in previous membership selections. We include 
the classical non members a-Cru and /3-Cru as new members as well as the pre-main- 
sequence stars HIP 79080 and 79081. We also show that the association is well mixed 
over distances of 8 degrees on the sky, and hence no determination can be made as to 
the formation process of the entire association. 

Key words: open cluster and associations: individual: Sco-Cen, stars: early-type, 
methods: statistical, stars: kinematics and dynamics, 



1 INTRODUCTION 

The Scorpius-Centaurus-Lupus-Crux OB Association (Sco 
OB2, Sco-Cen) is the nearest location with recent mas- 
sive star formation. The association was first identified by 



Kapteyn (19141 during an investigation of the parallaxes 



of 319 bright OB stars in the region of sky occupied by 
Sco-Cen. Following this, other kinematic studies confirmed 



that Sco-Cen is indeed a moving gro up (|Plaskett (19281, 



Blaauw| (|1946[) , |Bertiau| ( |1958[ ) , |Petrie| ( |1962[ ) , |Jones| ( |1971[ ) , 
de Zeeuw et al.| ( |1999| )). Since its discovery, Sco-Cen has 



been classically divided into three distinct sub-groups (see 
Figure [T]), Upper-Scorpius (US), Upper-Centaurus-Lupus 
(UCL), and Lower-Centaurus-Crux (LCC) ( |Blaauw|[T9i6| ) , 
with mean parallaxes of 6.9, 7.1 and 8.5 milli-arcseconds re- 
spectively, or distances of 145, 143 and 118 pc (de Zeeuw 
|et al.|[T999[ ). UCL and LCC have little interstellar material 
associated with them, whereas filamentary material can be 
observed towards US which is connected to the Ophiuchus 



cloud complex, a region of ongoing star formation ( de Geus 



1992a I. Photometry has demonstrated that the Ophiuchus 



cloud complex is on the near side of US at approximately 
125 pc, and isochrone fitting gives ages for the sub-groups 
as 5MyrforUS 16Myr for UCL and 16Myr for LCC (|3e 
Geus et al.|1989| ). 

US has received significantly more attention than the 
other two subgroups of Sco-Cen primarily due to its rela- 
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lively compact size. The age spread of US members is tightly 
bunched around 5Myr, which has led to the conclusion that 
star formation in US was an externally triggered process 
( Preibisch fc Mamajek|2008| ). The trig gering effect has been 
identified as the Shockwave created by several supernova ex- 
plosions in UCL which occurred approximately 12 Myr ago. 
Observations of the kinematics of the large HI loops sur- 
rounding Sco-Cen suggest that the Shockwaves passed US 
approximately 5 Myr ago, which agrees with the stellar age 
of the subgroup ( |de Geus|1992bl ). 

UCL and LCC have been subject to less study than 
US. This is primarily due to their relative lack of concen- 
tration on the sky and the closer overlap with the Galactic 
plane, which makes separation of members and field stars 
more difficult. Also, large portions of UCL and LCC are not 
observable from the northern hemisphere (where many early 
investigations were conducted), and early investigation into 
pre-main sequence stars focused on dark and reflection neb- 
ulae, which are less frequent in UCL and LCC ( |Preibisch| 
fc Mamajek||2008 1. Only in the last two decades have true 



high quaUty measurements of star motions on the sky (HIP 
PARCOS) become available, which have allowed separation 
of UCL and LCC members from the field \de Zeeuw et al 



1999 1. The formation processes which shaped UCL and LCC 



are as yet unknown, but are expected to be significantly 
more complicated than that of US. 



Elmegreen fc Lada ( 1977 1 propose that OB associations 



are formed through successive bursts caused by Shockwaves 
and ionisation fronts in molecular clouds. OB stars drive ion- 
isation fronts into neutral material in their vicinity. In dense 
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Figure 1. Positions of Sco OB2 on the sky with representative 
HIPPARCOS proper motion vectors. The members shown in this 
plot are taken from the |de Zeeuw et ah| ( |1999] l membership deter- 
mination, and show the windows (blue) used for each sub-group 
membership search. The familiar "Southern Cross" and the con- 
stellation Scorpius arc plotted in red, and the Coalsack nebula is 
shaded in blue. 



material the ionisation front is preceded by a shock front as 
it moves into a nearby molecular cloud. |Elmegreen fc Lada| 
1 1977 1 suggest that the gas accumulated between the two 



fronts becomes gravitationally unstable and can collapse to 
form stars. These new stars will then drive the process fur- 
ther into the remaining molecular cloud. This model agrees 
with the observations that OB associations tend to display 
a series of distinct subgroups ( Ambartsumian 19551. It is 



possible that Sco-Cen formed in this way. 

The latest high-mass membership study of Sco-Cen was 



carried out by de Zeeuw et al. ( 1999 1 who identified a total 



of 521 members. The search combined two different selec- 
tion methods, the convergent-point and "Spaghetti" meth- 



ods (de Bruijne 1999 Hoogerwerf & Aguilar 19991, which 



use proper motion and parallax to differentiate between 
members and non-members in distinct regions of sky corre- 
sponding to each Sco-Cen subgroup. The convergent-point 
method has long been a means by which to identify moving 
group members ( |Jones|1971[ ) , though it does have some clear 
weaknesses. The method is one-dimensional, with the selec- 
tion of member stars being based only on proper-motion 
in a specific direction. Additionally, the convergent-point 
method has a bias towards selecting more distant stars and 
stars with small proper motion. The "Spaghetti" method is 
stronger than the convergent-point method in that it takes 
into account both proper motion and parallax, though it 
does not constrain the radial velocity of members. This in- 
dicates that the development of a selection method which 
makes use of radial velocity is a clear avenue for improve- 
ment. 

The Sco 0B2 association is a useful and readily avail- 
able astrophysical laboratory. It is the largest group of newly 
formed stars in close proximity to the sun and provides 
an ideal testing ground for new star formation and evolu- 
tion models. Despite, this Sco-Cen is still relatively poorly 



studied when compared to other associations such as the 
Pleiedes, which is only a fraction of the size of Sco-Cen. 



MULTI-DIMENSIONAL BAYESIAN 
MEMBERSHIP SELECTION 



Since the last membership selection of |de Zeeuw et ah] 
( |1999[ ), there has been a vast improvement in the avail- 
able astrometric data for stars in the Sco 0B2 field of the 
sky. The HIPPARCOS ( [Perryman et al.|1997[ ) mission data 
has been reduced in a more accurate way by |van Leeuwen] 
(20071, and catalogues of radial velocities for a large num- 



ber of the Sco-Cen potential members have become available 
( jKharchenko et al.|2007[ ). These improvements allow one to 
develop a membership criterion based not only on more ac- 
curate proper motion but also parallax, radial velocity and 
galactic latitude. 

In the context of the membership analysis described 
in this paper, the HIPPARCOS data used was that from 
van Leeuwen (20071 of stars with B — V ^ 0.6, within 



the area of sky bounded by (285 ^ I ^ 360) degrees and 
(—10 ^ 6 5S 60). These are the bright, higher-mass objects in 
the area of sky occupied by the Sco OB2 association. The ra- 



dial velocities were taken from the Kharchenko et al. ( 2007 1 



catalogue (CRVAD-2) which is a cross-linking of the All-Sky 
Compiled Catalogue of 2.5 million stars ( Kharchenko|20(jT I , 



the General Catalogue of Radial Velocities ( Barbier-Brossat 



& Figon 2000 1 and a number of other, smaller radial ve- 
locity catalogues ( [Kharchenko et al.||2007] ). Measurements 



from the RAdial Velocity Experiment (RAVE) (Steinmetz 



et al. 2006 1 are not included in the CRVAD-2 catalogue. 



and as such, CRVAD-2 can be considered an up-to-date cat- 
alogue of stellar radial velocities in the pre-RAVE era. Given 
that RAVE observations of the Sco OB2 area of sky are not 
available, CRVAD-2 is the most complete source of radial 
velocities available for this work. Cross-checking was done 
between the HIPPARCOS and CRVAD-2 catalogues to find 
radial velocities for as many of the HIPPARCOS objects as 
possible. This resulted in 35% of objects having a known 
radial velocity, with typical uncertainties of 1 to 5 km s^^. 
Objects without a radial velocity were assigned a random 
radial velocity with an extremely large uncertainty, i.e. the 
radial velocity is treated as unconstrained in the Bayesian 
membership analysis described below. 

Our membership selection uses the six physical param- 
eters of each candidate star to decide upon a membership 
probability. These parameters are the three Galactic veloc- 
ity components, and the three position parameters, distance. 
Galactic longitude and Galactic latitude. In our new mem- 
bership selection, a member must have velocity. Galactic 
latitude and distance consistent with the five parameters 
derived from models described below. The five parameters 
derived from the models {U, V, W, r and b) will all be func- 
tions of Galactic longitude. 

The Galactic longitude constraint on membership is a 
simple cut at each end of the association. The Galactic lon- 
gitude of members is cut at / = 360 degrees at the US end 
of the association and a,t I — 285 degrees at the LCC end. 
At the US end, active star formation is in play so it is a 
logical end-point for the association. At the LCC end, mem- 
bership probabilities are typically lower because of the con- 
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Figure 2. The Galactic velocity components, distance and Galactic latitude of the |de Zeeuw et al.| ( [T999^ members. The linear trends 
with Galactic longitude are shown on each graph. Figure [2e] shows both the linear trend in Galactic latitude and the linear trend in the 
spread of Galactic latitude {cT},) as the dotted lines above and below the solid lino. The origin in Figure [2d| is the position of the sun, and 
stars with a greater distance have larger distance errors. 



fusion with the galactic plane, and the lack of very bright 
high-mass stars at lower galactic longitudes means there is 
no clear reason to extend our search further. There may of 
course still be associated recent star formation beyond LCC, 
but examining the possibility of this is beyond the scope of 
this study. In the following subsections we will describe our 



model for the association, the computational framework and 
the Bayesian membership selection. 
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Slope 


Intercept (1 = 0) 


u 


-0.13 


50.9 


V 


-0.009 


-20.8 


w 


-0.019 


-12.3 


b 


0.21 


-55.0 


Ob 


-0.01 


45.6 




A 


B 


r 


0.0059 


-0.0072 



Table 1. The linear model parameters derived for the |de Zeeuw| 
|et al.|p999| l members which were used. The galactic velocity (U , 
V , W) parameters are in units of km s~^, the galactic latitude 
(6, crj) parameters are in degrees and the distance (r) parameters 
are in pc~^ . 



2.1 Linear Models 

It is important to first note that the Galactic velocities of 



tlie de Zeeuw et al. ( 1999 1 members sliow a clear linear trend 
with Galactic longitude, with the trend being most pro- 
nounced in the component of the velocity pointing out of the 



Galactic Centre (U) (see Fig. 2a I. This trend is also present 
in the distance of group members, Galactic latitude and the 



spread of Galactic latitude of de Zeeuw et al. ( 1999 1 mem- 
bers. Note that the linear trend in the spread of latitudes 
was derived by assuming that the window sizes defining the 
subgroups (see boxes in Figure[l]) were a two-sigma coverage 
of the Sco-Cen association. This produced a standard devi- 
ation of Galactic latitude for each subgroup. A linear trend 
against Galactic longitude was then fitted to these three 
values. We denote the Galactic latitude standard deviation 
trend as (Ti,. 

We characterise the distance of the group members by 
a line in 3D space, which removes the anthropocentric bias 
associated with a simpler model where parallax is a function 
of galactic longitude . The line fitted to the distance of the 



de Zeeuw et al. ( 1999 1 group members is of the form; 



r = {A cos I + B sin I) 



(1) 



where A and B are the fitted parameters and I is the Galac- 
tic longitude. Figure [2d] displays this fitted line and the |de| 
Zeeuw et al. ( 1999 1 members. Note that the points with 
larger distances also have proportionally larger errors, and 
so these outlier points did not contribute to the fit in any 
significant way. 

The trends span all three subgroups (US, UCL and 
LCC) suggesting that perhaps Sco OB2 can be modelled 
as a continuum rather than three separate subgroups, with 
the galactic velocity components, distance and Galactic lat- 
itude following a linear trend. The linear trend fitted to the 



de Zeeuw et al. ( 1999 \ results will later be used in the cal- 
culation of membership probabilities. 

The IDL software package mpf it was used to obtain the 
parameters of the linear trends for the five parameter space 
variables. The fitted parameters can be seen in Table [l] Our 
choice of fitting a linear trend was based on there being no 
clear observational distinction between the subgroups, and 
is not a claim of continuous star formation. 



2.2 The Convergent-Point and Expected Values 

2.2.1 The Convergent-Pomt 

The membership selection method we have developed makes 
use of the concept of the convergent-point which was first 



introduced by Jones ( 1971 \ and modernised by de Bruijne 
( 1999 1. The convergent-point is the point on the sky at which 
a group of stars with a common Galactic velocity will ap- 
pear to converge to and has formed the basis of one of the 
search methods used to populate the current membership 
of Sco-Cen (de Zeeuw et al. 19991. The underlying princi- 



ple of our new membership selection and its relation to the 
convergent-point is as follows: For each candidate member, 
expected Galactic velocity components {Ug{l), Vg{l), Wg{l)), 
distance (rg) and Galactic latitude (bg) can be calculated 
from the linear trends as functions of the Galactic longitude 
(l) of the candidate. From the expected Galactic velocities, 
a convergent-point is calculated {lcp,bcp) and based on this 
convergent-point, expected proper motions and radial veloc- 
ity can be found (described below). Hence, for a given star 
to be a considered a member, it is required to have proper 
motions, radial velocity and distance which are consistent 
with the expected values calculated from the linear models 
in Galactic longitude. 

The significance of the convergent-point to the member- 
ship selection developed in this paper relates to the velocity 
coordinate system of the data. The velocities provided by 
the HIPPARCOS and CRVAD-2 catalogues are two proper 
motions (equatorial coordinates) and one radial velocity, 
while the linear models provide Galactic velocity and paral- 
lax as a function of Galactic longitude. The concept of the 
convergent-point allows a direct and easily computed con- 
version from the model Galactic velocities to the velocity co- 
ordinate system of the available data. The convergent-point 
can be calculated as follows; 



tan(Zc 



tan(6cp) = 



(2) 



where Icp and hep are the coordinates of the convergent-point 
on the sky and Ug{l), Vg{l), and Wg{l) are the three Galactic 
velocity components derived from the linear trends and are 
functions of Galactic longitude. Thus the convergent-point 
is also a function of Galactic longitude. 

Once the convergent-point for a given candidate mem- 
ber is found, the proper motions provided by the HIPPAR- 
COS catalogue are rotated into a new coordinate system 
with proper motion along the great circle joining the star 
and the convergent-point and proper motion along the great 
circle perpendicular to this. The angle of rotation (7) is 
shown in Fig. [s] (in Galactic coordinates) and is calculated 
using trigonometric identities. Note that the rotation angle 
can also be calculated in Equatorial coordinates and will be 
referred to as 9. The Equatorial rotation angle 6 was used 
to rotate the HIPPARCOS proper motions into the new co- 
ordinate system described above. The rotation is given by 
the following equation; 



sin 9 

- cos 6 



cos ( 
sin I 



(3) 
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Figure 3. The triangle made from the three great-circles joining 
the star, convergent-point and Galactic north pole. The angles 7 
and A are the rotation angle and angular distance between the star 
and the convergent-point respectively. Note that coordinates used 
here are Galactic coordinates, and also that an analogous diagram 
can be drawn for equatorial coordinates with the rotation angle 
9 and the same angular distance A. 



where fj,c, and are the two HIPPARCOS proper mo- 
tions, /ix is the proper motion perpendicular to the direction 
of the convergent-point and is the proper motion towards 
the convergent-point. 

The use of the convergent-point in our method allows 
the parameters of the models to be easily converted to the 
framework of the observational data, rather than manipu- 
lating the data into a form comparable to the models. The 
method also allows uncertainties to be calculated via a sim- 
ple rotational transformation of the HIPPARCOS covari- 
ance matrix rather than complicated conversions to Galac- 
tic velocity coordinates. It is also important to stress that 
the convergent-point is a function of Galactic longitude (al- 
though indirectly) and hence each star will have a different 
convergent-point. Hence, we are not using the convergent- 
point in the traditional sense as the point at which asso- 
ciation members converge on the sky, but as a means of 
constructing a new coordinates system which allows us to 
convert from Galactic velocities to proper motions and ra- 
dial velocity. Furthermore, this convergent-point based coor- 
dinate system allows a simple membership criterion for /xx: 
this quantity should be small for a member. This criterion is 
the basis for the classical convergent-point selection method 



used by Jones ( 1971 1 



dial velocity require conversion from the linear models for 
Galactic velocity components {Ug{l), Vg{l), Wg{l)). 

To convert the expected Galactic velocity from the lin- 
ear models, unit vectors in the direction of and the ra- 
dial velocity of each star are needed. These will be denoted 
p and f respectively and can be calculated using spherical 
trigonometry; 



p = cos 7b -I- sin7l, 
(cos I sin b, sin I sin fe, cos /), 



(4) 



where 1 and b are the Galactic longitude (l) and Galac- 
tic latitude (b) unit vectors. The expected Galactic velocity 
{ug{l)) is then projected in the direction of these two unit 
vectors yielding the expected velocity along the great circle 
joining the convergent-point and the star, and the expected 
radial velocity, denoted uii and Vr ■ 



Urg = f ■ Ug{l). 



(5) 



In order to compare the expected parallel velocity (i/y ) 
to the proper-motion in the direction of the convergent-point 
for a given star, it is required that be converted into 
proper-motion units. This was done by multiplying by the 
parallax(n) and a M = 4.74 kmyrs"^ which is the ratio 
of the number of seconds in a Julian year and one astro- 
nomical unit in kilometres, yielding the expected parallel 
proper-motion for the given star. This allows ideal member- 
ship to be described by Equation Q . Note that the expected 
quantities will vary depending on the position of the star in 
question, and so will be different for each possible member. 



MX 







(6) 



= HMi/||g(/cp,bcp) 

r = Tg{l) 
b = bg{l) 

The above five conditions describe the expected param- 
eters of an ideal member star at a given Galactic longitude. 
It is expected that there is a velocity dispersion around the 
expected Galactic velocity of each star due to gravitational 
interaction of Sco-Cen members. The magnitude of the ve- 
locity dispersion used in this selection method was 3 km s~^ 
( Kraus fc Hillenbrand|2008 |. 



2.2.2 Expected Values 

The basis of this new membership selection is that for a 
candidate star to be a member it must have velocity, dis- 
tance and position which are consistent with some expected 
values. The expected values of the five parameters are de- 
rived from the linear model described above, and vary for 
each candidate star with Galactic longitude. The expected 
Galactic latitude and parallax are taken directly from the 
linear models, while the expected proper motions and ra- 



2.3 Preliminary Exclusions 

Using the expected velocities expressed above and the er- 
rors associated with the Hipparcos and radial velocity data, 
a membership probability can be calculated for each star, 
using Bayes' factors and Monte-Carlo integration. Before 
the probabilities are calculated, a simple and preliminary 
yes/no selection was made to initially reject stars which 
would have low membership probabilities. This was done by 
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finding the squared, inverse- variance weighted difTerence be- 
tween each velocity component and its respective expected 
velocity (Equation ([7|), and using a three degree-of- freedom 
probability density function to calculate the probability, 
?7, of the sum being greater than its calculated value (Equa- 
tion (|8|). Stars for which a radial velocity was not available 
in the Karchenko catalogue were treated with a two degree- 
of-freedom probability density function, with the radial 
velocity term R = (see Equation ([7|). Ideally, model like- 
lihood ratios would be calculated for each star in the field, 
however, stars with a large probability rj will have a small 
Bayes' factor and will be excluded by the Bayesian process 
anyway. 



two models, in the case of this membership selection the 
models are that a given star is a member of Sco-Cen, denoted 
Mg, and that a star is a member of the field and not a 
member of Sco-Cen, denoted Mf. The probability of a model 
given the data is R/{R + 1). The model likelihood ratio 
(R) is given by Equation ^ ( |Sivia fc Skilling||2007| ), where 
D denotes the data, the vertical bar denotes "given" and P 
denotes probability. 



R = 



P{M,\D) 
P{Mf\D) 



(9) 



Using Bayes' Theorem (Sivia & Skilling 20071 the model 



likelihood ratio can be expressed as: 



T : 



P = 



(M||-Af,.|l^n)^ 



(7) 



7? = 



t2 a„2 I 



V = 



2iT{ 



L 



x^e 2 da;. 



(8) 



T+P+R 



The cut off value used for 77 for stars which progressed 
to the Bayesian analysis was determined by testing whether 
setting T) lower would change the number of stars chosen by 
the Bayesian analysis at a 50% threshold. A value of 0.8 was 
used for r) , which was well above the condition described, 
but still able to significantly reduce the time required for 
computation. Within Sco-Cen, there is a common velocity 
among members and a dispersion away from the common 
velocity caused by gravitational interaction. In Equation 
([7| we have included internal velocity dispersion Gint as a 
source of deviation from the expected velocities. Conversion 
to proper motion units is required for /xy and using the 
conversion factor M = 4.74 km yrs~^ mentioned above and 



the mean parallax of the de Zeeuw et al. ( 1999 1 members 



The converted velocity dispersion is denoted Oint, ■ 

This preliminary selection did not produce membership 
probabilities, but simply functioned to filter out stars which 
could not be members based on their velocity characteris- 
tics. This is why we do not present a comprehensive analysis 
of the stars selected, though a brief account is as follows. Ap- 
proximately four times the number of stars in the |de Zeeuw| 
et al. ( 1999 1 membership were included by this preliminary 



selection (2051 stars), which is expected since it was de- 
signed to let through all but the most clear non-members 



(5473 excluded). Note that 15 |de Zeeuw et al.| ( |1999[ ) mem- 
bers were rejected, almost all of which were due to a ra- 
dial velocity inconsistent with membership, again, this is 
expected as radial velocity was not used by |de Zeeuw et aT] 
( |1999[ ). A further discussion and additional details on these 
stars will be given in the following sections. 



2.4 Bayesian Membership Probabilities 

To find a membership probability for each star in the field of 
Sco OB2, Bayesian model testing was used. Bayes Factors 
can be used to calculate the likelihood ratio (R) between 



R ■- 



P{M,) P{D\M,) 



K, 



(10) 



PiMf) P{D\Mf) P[Mj) 
where K is the Bayes' factor, and can be written as 



P[D\M,) ^ SP{4>,\M,)P{D\ct>^)Act,^ 
P{D\Mf) S P{4>f\Mj)P[D\4>f)d4>f' ^ ' 

where <f)g and are a set of parameters {[/, V, W, r, b}. 

The components which make up the five dimensional 
integrals in equation ( |11[ ) can be separated into three un- 
coupled integrals; 



PiD\M,) = P{b,\M,)P{D\b,)db, 



r 

•J o< 



P{r,\M,] 



P{D\rg)Ar, 



(12) 



/P({U}lAfg)P(D|{U})d{U} 



where {U} is the set of three group Galactic veloc- 
ity components {Ug,Vg,Wg}. These cannot be separated 
because the data do not constrain these components indi- 
vidually. Note that a completely analogous equation can 
be written for the field integral (denominator of equation 
(111). Firstly, in the Galactic latitude (b) integral, the term 
P[D\bg) is the delta function 5{b — bg) because the uncer- 
tainties in the HIPPARCOS positions are trivially small. 
Hence the Galactic latitude integral is analytically computed 
to be simply P{b\Mg), which is a Gaussian with mean and 
standard deviation taken from the linear models, and trun- 
cated at the edge of our search area (— 10 < 6 < 60). Note 
that since the mean and standard deviation are function 
of Galactic latitude, the Gaussian will be different for each 
star. Similarly, for the field, P{D\Mf) is a Gaussian with a 
mean of zero and standard deviation taken from the spread 
of Galactic latitudes in the HIPPARCOS catalogue. Again, 
this was truncated at the edges of our search area. 

A more complicated issue is that of the distance (r) in- 
tegrals. It is important to note that there is a population 
effect in the HIPPARCOS catalogue towards more distant 
stars. This is expected as there are intrinsically more stars 
per unit distance from the sun. This effect will thus be in- 
trinsic to both group and field distance distributions. Hence, 
the field distribution will be defined by this population ef- 
fect, and the group distribution will be treated as a Gaussian 
superimposed on the population effect. The probability dis- 
tribution of the population effect was determined by binning 
all the HIPPARCOS stars in our search area into 5pc bins 
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and calculating the fraction of stars per bin. The population 
distribution will be denoted B{r). The terms in the distance 
integral components for the field are then given by the fol- 
lowing; 



P{(j)f\Mf) oc expl 



I (v-Vfr 

+ 2al 



(18) 



P(rf\Mf) = B{r) 



(13) 



(w-WfY 



)) 



P{D\rf) = 



exp I 



)■ 



(14) 



The integral of these two components is computed nu- 
merically by first defining the 5 pc distance bins and calcu- 
lating the fractional number of stars per bin. This produces 
a set of r/ values and P(r/|M/), their weightings (one for 
each bin). For each of these r/ values P{D\rf) is calculated, 
multiplied by it's corresponding weighting and then totalled 
to yield the integral. The corresponding group distance dis- 
tributions are defined in a similar way; 



P(r,\M,) = P(r,\Mf)Ce^p{-{ 



P{D\r, 



exp I 



jr.- foil))' 



2a2 



)), (15) 



(16) 



where fg{l) is the distance taken from the linear trend 
with Galactic longitude, and C is a normalisation constant 
which is computed numerically, af is the standard deviation 



of the de Zeeuw et al. ( 1999 1 members used to fit the linear 



distance trend. The assumption that the |de Zeeuw et al.] 
1 1999 1 members describe the depth of the association well is 



valid because distance was not constrained in the de Zeeuw 



et al. (19991 membership selection. Note that, as described 



above P{rg\Mg) is a Gaussian in distance superimposed on 
the population distribution, and the factorP(D| rg) will be 
different for each star since fg{l) is a function of Galactic 
longitude I. The group distance integral is computed in the 
same way as the field integral. 

Having separated the Galactic longitude and distance 



components from the integral in Equation (111, the remain- 



ing integral is the contributions from the Galactic velocity 
components. The factors involved in this integral are de- 
scribed below. 

P{(f)g\Mg) (The probability of the parameters given the 
model that a star is in the group) is modelled as the prod- 
uct of three Gaussian distributions and is shown in Equation 
( |17[ ). The means of the Gaussians are the expected values 
of the parameters derived from the linear trends in Galactic 
longitude {Ug{l),Vg{l) and Wg{l)). The spread of these pa- 
rameters around the expected values is only influenced by 
the internal velocity dispersion of the association, thus the 
Gaussians have a standard deviation equal to atnt- 



P{(i>g\Mg) oc expl 



The distribution for the field, P{tpf\Mf) (Equation 
( |18[ )), is defined in an analogous way. However, the Galactic 
velocities are taken from the Galactic Thin Disk with values 
of ([//, Vf,Wf) = (9.0, -6.9, -7.0) km s"^ and standard de- 
viations of auj = 19.8, avf = 12.8, and awf = 8.0 km s~^ 



for A-type stars ( [Robin et aL]|2003 |. 

The probability of the data given the parameters from 
the group, P{D\cf>g), is also defined as the product of three 
Gaussian distributions, though rather than Galactic veloc- 
ity, the new coordinates calculated from the convergent- 
point using the projections in equation ([5| were used, i.e. 
, fj,± and Vr- 

The final remaining probability distribution to define 
is P(D|(/)y^), the probability of the data given the model 
that the star is in the field. This was also modelled as three 
Gaussians with expected values of /ix, /iy and Vr calculated 
in the same way as was done for the group. Note that this 
involves finding a convergent-point for various points in the 
parameter space. One may ask why the field stars have a 
convergent-point, and the answer is that the apparent mo- 
tion of the field stars is dominated by the motion of the sun, 
and the rotation of the galactic disk, and so a convergent- 
point is in fact present for a star at rest with respect to the 
Galactic Disk. However, in this case the convergent point 
is really an abstract coordinate system used to find and 
. Hence expected values are calculated for each star as de- 
scribed in the previous section. These two probability func- 



tions are shown below in Equations ( 19 1 and ( 20 1 



P(Djc/)„) oc exp 



+ 



MX 

(M||-nMi.|, (i) f 



2(<^J|l+(Mi.||g(i)<Tn)2) 



(19) 



P{D\c^f) (X exp(-(^ 

*'x 

(M||-AAi||,n)2 



2('^J|l+('4M||^<Tn)2) 



(20) 



+ ^)) 



-I ^ 



{w-Wg(i)r 



)) 



(17) 



The above was implemented using a Monte-Carlo 
scheme to compute the integrals in equation ( |11[ ). 10''' sam- 
ples were used for the group integral and 10^ were used for 
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the field to adequately sample the larger parameter space of 
the field distribution. A brief overview of the process is as 
follows: From the model distributions {Mg and M/), a set of 
randomly sampled parameters are taken to form (p^ and 4>p 
Mg is in fact a normal distribution in each parameter (U, V, 
W, r, b), with means taken from the linear trends in Galac- 
tic longitude. Hence, cj)^ — {U,V,W,r,b} is a function of 
Galactic longitude and will be different for each star. From 
each parameter sample, a convergent-point is calculated and 
the expected values of jj,±, /iy and Vr are found with respect 
to each convergent point. Once this is done, P{D\(f)g) and 
P{D\(f>j:) can be found and the integral can be calculated. 
Note that the integral is over the parameters U, V, and W 
and not over Galactic longitude, however, these parameters 
are functions of Galactic longitude for the group integral. 



COS proper motions. Of the latter 53, only 9 were excluded 
based on perpendicular proper motion {n±), which was a 



The prior used in the Bayesian selection method 



' P(Mf ) ' 



was generated in an iterative manner. The first prior was 
the ratio of the number of member stars in the search area 
to the number of non-member stars in the lde Zeeuw et al.l 
( 1999[ ) catalogue, with B-V colour less than 0.6 magnitudes. 
An A-type star with B-V of 0.6 magnitudes at the mean 
distance of the Sco-Cen association is expected to appear to 
be 9"" magnitude in the visual. Beyond this magnitude the 
completeness of the HIPPARCOS data greatly declines and 
so selection of fainter members is not feasible. Subsequent 
priors were taken as the ratio of the total of the probabilities 
of membership and the total number of stars in the field. 
This was repeated until the prior converged (ten iterations) . 
The prior iteration is described in the following equation, 
with pn being the membership probabilities, Pk being the 
fc*'' prior and Nf being the number of stars in the field. 



PiMg) 



(21) 



Once the prior converged, stars with membership prob- 
abilities of 50% or greater were denoted members of the 
Sco-Cen association, which has been the historical conven- 
tion. 



major basis of the de Zeeuw et al. (19991 selection. The ex- 
clusion of previous members based on HIPPARCOS proper 
motions, especially the exclusions based on px, emphasises 
the significant effect of the newly reduced data in the mem- 
bership selection. Nine of the rejected previous members 
(HIP 62002, 64661, 68534, 72099, 74449, 77315, 79610, 80208 
and 81316) had membership probabilities greater the 40%, 
but did not have accurate or available radial velocity mea- 
surements. It is likely that these stars can indeed be shown 
to be members of Sco-Cen with accurate radial velocity mea- 
surements. 

Importantly, the membership selection has established 
the classical non-members a-Cru and /3-Cru as Sco-Cen stars 
with membership probabilities of 66 and 73 % respectively. 
This is the first time these stars have been included as mem- 
bers of Sco-Cen. It is important that they are included as 
they are bright B-type stars with the correct position and 
parallax and age to have formed as Sco-Cen members. 

Also, we select the classical members r;-Cen and fi^-Sco 
which were not included in the last membership survey, and 
we confirm the membership of HIP 59449, as suggested but 
unconfirmed as a Sco-Cen member ( de Zeeuw et al.|1999 l. 

De Zeeuw et al. (1999) suggested, but was unable to 
confirm that the pre-main-sequence stars HIP 79080 and 
HIP 79081 were Sco-Cen members. These stars have now 
been selected as members with probabilities of 51 and 63 %. 
HIP 78094, another pre-main-sequence star, has consistent 
proper motions for membership but cannot be confirmed 
as a member due to a poor parallax measurement. We also 
confirm that the classical member 5-Sco has proper motions 
which are inconsistent with group membership. This star 
was also excluded from Sco-Cen by de Zeeuw et al. (1999), 
with the explanation that its binarity, with period of ~ 20 
years could be responsible for producing a proper motion 
inconsistent with the group. 



3 THE NEW HIPPARCOS MEMBERSHIP OF 
SCO-CEN OB2 

The membership selection method developed in this paper 
was used to produce an improved membership list for the 
Sco-Cen OB association. A full list of member stars and 
their probabilities can be seen in Appendix |A] Figure |4] dis- 
plays the selected members on the sky with proper motion 
vectors. Members of Sco-Cen have now been identified in a 
much larger area on the sky, with new members at signifi- 
cant distances from the central clustering of the association. 
While there is clear distinction between Upper-Scorpius and 
the rest of the association, UCL and LCC now show consis- 
tent density and even clumping of members, including in 
regions that were previously considered sparser boundaries 
between subgroups. 

The selection has produced 436 Sco-Cen members, 348 
of which were previous members ( de Zeeuw et al.|1999 \ with 
the rest being new members. 93 de Zeeuw et al. ( 1999 1 mem- 



4 DISCUSSION 

The new membership described here does not require the 
segmentation of Sco-Cen into the three classical sub-groups 
in order for members to be chosen. Instead Sco-Cen has been 
treated as a continuum using a linear model to select mem- 
bers and reject non-members. The use of a linear model has 
removed the need for hard-edged boundaries around selec- 



bers have been rejected as members of Sco-Cen, 11 on the ba- 
sis of parallax, 29 from radial velocity and 53 from HIPPAR- 



tion areas, as seen in the previous membership (de Zeeuw 
et al.|1999 l, and has allowed the area in which members can 
be found to be slightly expanded. Furthermore, the linear 
trend in the component of the Galactic velocity pointing out 
of the Galactic Centre (U) demonstrates that the association 
is indeed expanding. This is expected since OB associations 
are considered to be loosely bound and is consistent with 
the idea of a common origin for the entire association. 

Fig ure [5| displays phot ometry taken from the Tycho cat- 
alogue ( [Perryman fc ESA|p97l ) in a colour-magnitude dia- 
gram of the new Sco-Cen members. The Tycho photometry 
has been converted to standard UBV using the table found 
in |Bessell| (|2000[) . We have plotted three Padova isochrones 



( Girardi et al.||2002| at 4, 10 and 20Myrs ages which have 
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Figure 4. Position and proper motion vectors of tlie new Sco-Cen members. Note the significantly expanded region of sky within which 
members have been selected. Diamond shaped points are stars which were included in the last membership | |de Zeeuw et ah|1999[ l and 
triangular points are new members which were not present in the last membership compilation. Darker blue indicates higher membership 
probability 
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Figure 5. Colour-magnitude diagram of the Sco-Cen members selected by the new method with Padova isochroncs for 4 (solid), 10 
(dashed) and 20 (dotted-dashed) Myrs. Note the two isolated stars, these are shifted due to differential reddening caused by the p-Oph 
cloud complex which impinges on parts of Upper Scorpius. 
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Table 2. List of stars separated from the main sequence in Fig. 



HIP 


I 


b 


HIPPARCOS Spectral Type 


78820 


353.2 


23.6 


B0.5V 


78933 


352.7 


22.8 


BIV 


79374 


354.6 


22.7 


B2IV 


80112 


351.3 


17.0 


BlIII 


80569 


357.9 


20.7 


B2Vne 



Distonce (deg) 



Figure 6. A logaritlimic histogram of the two point correlation 
function applied to the new Sco-Cen membership list. 



been reddened to match the Tycho photometry using the 
following equation ( |Aumer fc Binney|2009| . 



E{B - V) = 





0.47(d- 



70pc)/lkpc 



d < 70pc 
d > 70pc 



(22) 



The value of E{B — V) was calculated to be 0.027 mag- 
nitudes and the isochrones were reddened by this amount. 
From the HR diagram, it can be seen that the youngest 
stars in Sco-Cen are approximately 6 Myr old, while the 
older stars appear to be up to 20 Myr old. This is in agree- 
ment with historical estimates of the age of the association 
| Blaauwl|1964 1. Also, Upper Scorpius, the youngest section 
of Sco-Cen ,is expected to be have formed approximately 
5 Myr ago ( Preibisch et al.|2002 1 , which agrees with the age 
estimate of the youngest stars in the new membership. Note 
that the large spread of members around the isochrones in 
Fig. [5] is due to both the parallax measurement error, the 
effects of unresolved binarity on the photometry and the age 
spread. There are five stars in Figure [5] which appear to be 
above the main sequence, these stars are listed in Table [2] 
These five stars, based on their spectral types, are signifi- 
cantly more reddened than the rest of the new association 
members. This could be due to the impingement of the p- 
Oph clouds between the sun and Sco-Cen jde Gcus 1992a|. 

A two point correlation function was also applied to 
the new members of Sco-Cen for the purpose of examin- 
ing the effects of internal velocity dispersion on mixing the 
primordial fractal structure of the association. This meant 
totalling the number of other Sco-Cen stars within certain 



distance bins. The distance bins were: (0, 0.125), (0.125, 0.5), 
(0.5,1), (1,2), (2,4) and (4,8) degrees. Each bin was then 
area normalised. A histogram of the two point correlation 
function can be seen in figure [6] At small separation scales 



bound binary pairing produces a steep power law (Kraus 



& Hillenbrand 2008). Beyond separations of half a degree, 
random motions appear to blur out any primordial group 
structure and produce a shallow slope in the diagram. In 
this study it is seen that structure on the order of eight de- 
grees is blurred out by the group's velocity dispersion. This 
implies that hard-edged boundaries between the subgroups 
of Sco-Cen cannot be defined in a non-arbitrary way. Addi- 
tionally, this means that it cannot be determined if Sco-Cen 
formed via the shock front mechanism of [Elmegreen &: Lada| 
(119771, as evidence of separated structures has been effec- 



tively mixed out. 



4.1 Future Work 

The new membership of Sco-Cen presented here is for the 
bluer and higher mass population of Sco-Cen, with a colour 
B — V < 0.6. These are the brightest and most well known 
stars in Sco-Cen. It has been proposed that there are many 
thousands of lower mass stars which are members of Sco-Cen 
but which are yet to be identified. The membership selection 
method developed here can be directly applied to these lower 
mass stars once the required data becomes available. The 
only difficulty is deciding upon the use of Bayesian prior. In 
this study, we have used the |de Zeeuw et aL] ( |1999"l ) member- 
ship as a starting point in defining the prior which resulted 
in the new membership. There is currently no equivalent for 
the fainter, lower mass members. One possibility is to ex- 
trapolate the lower mass population based upon the higher 
mass population of the new membership in Appendix B, and 
use this extrapolation as a new prior. This would involve the 
use of accurate models to provide a colour-mass relation for 
young stars. Along a similar line, spectroscopic techniques 
can be used to identify lower mass Sco-Cen members. This 
has been done for a small region of Sco-Cen by IPreibisch] 
et al. (20021 using the multi-object spectrograph 2dF, but 



can be greatly extended with the use of an instrument such 
as 6dF at the UK Schmidt telescope which is better matched 
to the large area of sky covered by the association. As a final 
point on future membership work, the membership selection 
method developed here, with some adaptation, can be ap- 
plied to other OB associations and moving groups, as was 
done by |de Zeeuw et aL] ( |1999 | with their selection method. 

Age estimation for the various regions of Sco-Cen has 
traditionally been achieved through isochrone fitting. How- 
ever, due to parallax uncertainties and the difficulty in deter- 
mining the temperature of individual stars accurately, pre- 
vious age estimates have typically been based on a small 
number of stars ( [Preibisch et al.||2002[ ). Determination of 
dynamical masses for a number of Sco-Cen targets can al- 
low a more accurate age estimate, and more precise dating 
across the association. 
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APPENDIX A: MEMBERSHIP TABLE 

Table Al: Membership List and Probabilities 



HIPPARCOS Number Probability (%) 



50520 


62 


50847 


62 


52357 


82 


52742 


61 


53701 


68 


53771 


91 


54231 


88 


54767 


64 


55188 


78 


55334 


58 


55425 


61 


55899 


72 


56227 


58 


56354 


84 


56379 


83 


56963 


81 


57375 


51 


57710 


81 


57809 


82 


57851 


73 


57947 


50 


57950 


76 


58146 


77 


58167 


87 


58220 


88 


58283 


53 


58416 


87 


58452 


82 


58465 


88 


58528 


85 


58720 


59 


58859 


83 


58884 


78 


58899 


84 


58901 


51 


59119 


60 


59173 


66 


59257 


63 


59282 


88 


59397 


86 


59413 


86 


59449 


60 


59481 


88 


59502 


86 


59505 


81 


59603 


82 


59693 


74 


59716 


64 


59724 


59 


59747 


75 


59781 


82 


59898 


66 


59960 


89 


60009 


79 


60205 


62 


60245 


74 
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Table Al: Membership List and Probabilities 



HIPPARCOS Number Probability (%) 



60348 


78 


60379 


53 


60459 


88 


60513 


76 


60561 


73 


60577 


81 


60710 


71 


60718 


66 


60823 


79 


60851 


82 


60855 


70 


61049 


85 


61087 


88 


61257 


84 


61265 


67 


61426 


54 


61498 


81 


61585 


71 


61639 


84 


61684 


73 


61782 


84 


61796 


61 


62026 


77 


62032 


73 


62058 


80 


62134 


83 


62171 


78 


62179 


64 


62327 


84 


62427 


65 


62428 


81 


62431 


59 


62434 


73 


62657 


83 


62683 


62 


62786 


55 


63003 


68 


63005 


81 


63007 


83 


63041 


83 


63204 


87 


63210 


75 


63272 


84 


63439 


70 


63836 


85 


63839 


73 


63886 


75 


63945 


79 


63975 


57 


64004 


73 


64044 


78 


64053 


82 


64184 


81 


64264 


58 


64320 


70 


64322 


74 


64425 


54 


64515 


71 
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Table Al: Membership List and Probabilities 



HIPPARCOS Number Probability (%) 



64617 


51 


64837 


69 


64877 


61 


64892 


74 


64925 


80 


64995 


79 


65021 


59 


65089 


84 


65112 


81 


65136 


63 


65178 


61 


65219 


78 


65271 


62 


65348 


50 


65394 


67 


65426 


76 


65517 


69 


65617 


64 


65822 


81 


65875 


88 


65965 


70 


66068 


75 


66075 


56 


66447 


93 


66454 


89 


66566 


82 


66651 


75 


66722 


92 


66821 


83 


66908 


88 


67036 


81 


67068 


80 


67199 


86 


67230 


65 


67260 


61 


67428 


73 


67464 


57 


67472 


54 


67497 


92 


67669 


66 


67703 


70 


67919 


66 


67957 


89 


67970 


90 


67973 


75 


68080 


60 


68245 


76 


68282 


77 


68335 


87 


68413 


66 


68532 


89 


68722 


72 


68781 


89 


68862 


58 


68867 


56 


69011 


77 


69113 


80 


69291 


88 
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Table Al: Membership List and Probabilities 



HIPPARCOS Number Probability (%) 



69302 


82 


69475 


59 


69618 


71 


69720 


83 


69791 


60 


69995 


70 


70149 


90 


70300 


86 


70350 


70 


70441 


93 


70455 


51 


70483 


52 


70558 


88 


70626 


79 


70689 


90 


70697 


89 


70753 


80 


70809 


68 


70833 


84 


70904 


66 


70931 


72 


70998 


89 


71140 


67 


71271 


55 


71321 


92 


71352 


82 


71353 


90 


71453 


66 


71536 


84 


71708 


84 


71724 


83 


71860 


70 


71865 


93 


72149 


51 


72158 


59 


72216 


57 


72584 


78 


72627 


93 


72683 


74 


72800 


71 


72940 


91 


73145 


91 


73147 


78 


73150 


70 


73266 


92 


73334 


73 


73341 


65 


73535 


62 


73559 


61 


73624 


86 


73742 


58 


73807 


76 


73913 


88 


73937 


86 


73990 


81 


74066 


67 


74098 


67 


74100 


87 



A.C. Rizzuto et al. 

Table Al: Membership List and Probabilities 



HIPPARCOS Number Probability (%) 



74104 


74 


74468 


70 


74479 


78 


74499 


90 


74651 


62 


74865 


90 


74911 


62 


74950 


57 


74959 


84 


74985 


81 


75035 


61 


75056 


84 


75077 


78 


75151 


81 


75164 


78 


75210 


91 


75304 


90 


75367 


80 


75476 


82 


75480 


89 


75491 


70 


75509 


92 


75613 


52 


75647 


83 


75683 


83 


75802 


50 


75824 


69 


75891 


90 


75915 


92 


76001 


73 


76048 


83 


76063 


55 


76071 


70 


76084 


87 


76126 


58 


76143 


51 


76234 


73 


76297 


91 


76310 


69 


76371 


74 


76395 


78 


76501 


50 


76591 


64 


76600 


53 


76633 


89 


76875 


89 


76945 


85 


76997 


62 


77038 


88 


77086 


83 


77124 


56 


77144 


79 


77150 


90 


77286 


84 


77295 


74 


77317 


67 


77347 


56 


77388 


78 
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Table Al: Membership List and Probabilities 



HIPPARCOS Number Probability (%) 



77399 


59 


77432 


87 


77502 


62 


77520 


90 


77545 


84 


77562 


70 


77635 


89 


77815 


79 


77840 


51 


77858 


83 


77859 


75 


77900 


82 


77909 


71 


77911 


83 


77960 


86 


77968 


59 


78043 


83 


78099 


88 


78104 


82 


78150 


84 


78183 


88 


78196 


90 


78207 


86 


78233 


64 


78263 


59 


78324 


78 


78384 


84 


78416 


74 


78494 


54 


78530 


76 


78533 


88 


78549 


89 


78555 


89 


78581 


86 


78641 


91 


78655 


90 


78663 


88 


78702 


90 


78711 


63 


78754 


73 


78756 


75 


78809 


91 


78820 


50 


78821 


67 


78830 


63 


78847 


89 


78853 


55 


78877 


61 


78918 


72 


78933 


91 


78956 


93 


78963 


72 


78968 


89 


78996 


86 


79031 


88 


79044 


81 


79054 


91 


79078 


58 



A.C. Rizzuto et al. 

Table Al: Membership List and Probabilities 



HIPPARCOS Number Probability (%) 



79080 


51 


79081 


63 


79097 


66 


79098 


55 


79124 


82 


79142 


54 


79156 


88 


79229 


76 


79250 


92 


79288 


87 


79363 


71 


79366 


87 


79369 


82 


79374 


72 


79383 


74 


79392 


86 


79399 


82 


79400 


54 


79404 


82 


79410 


79 


79439 


63 


79476 


80 


79516 


77 


79530 


79 


79599 


85 


79622 


76 


79631 


86 


79673 


83 


79690 


70 


79710 


79 


79739 


65 


79742 


76 


79771 


89 


79785 


77 


79860 


72 


79878 


89 


79897 


90 


79910 


88 


79977 


90 


80024 


82 


80036 


68 


80059 


84 


80062 


57 


80088 


87 


80112 


56 


80126 


86 


80130 


89 


80238 


92 


80311 


86 


80324 


92 


80371 


85 


80390 


79 


80425 


66 


80461 


81 


80475 


52 


80493 


86 


80535 


90 


80569 


55 



Multi- Dimensional Bayesian Membership Analysis of the Sco 0B2 Moving Group. 



Table Al: Membership List and Probabilities 



HIPPARCOS Number Probability (%) 



80591 


84 


80711 


53 


80799 


91 


80815 


62 


80896 


86 


80911 


85 


81208 


58 


81266 


89 


81447 


59 


81455 


76 


81474 


91 


81477 


53 


81624 


79 


81751 


78 


81891 


58 


81914 


58 


81949 


74 


81972 


59 


82069 


51 


82218 


89 


82319 


51 


82397 


81 


82430 


78 


82514 


65 


82534 


89 


82545 


86 


82554 


71 


82560 


75 


82569 


63 


82714 


54 


83159 


62 


83508 


64 



